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Abstract. The universe was reionized by redshift z ss 6 by a small fraction of the 
baryons in the universe, which released energy following their condensation out of a 
cold, dark, and neutral IGM into the earliest galaxies. The theory of this reionization 
is a critical missing link in the theory of galaxy formation. Its numerous observable 
consequences include effects on the spectrum, anisotropy and polarization of the cosmic 
microwave background and signatures of high-redshift star and quasar formation. This 
energy release also created feedback on galaxy formation which left its imprint on the 
mass spectrum and internal characteristics of galaxies and on the gas between galaxies 
long after reionization was complete. Recent work suggests that the photoevaporation 
of dwarf galaxy minihalos may have consumed most of the photons required to reionize 
the currently-favored ACDM universe. We will review recent developments in our 
understanding of this process. 



I INTRODUCTION 

Observations of quasar absorption spectra indicate that the universe was reion- 
ized prior to redshift z = 5 (e.g. [68]). CMB anisotropy data on the first acoustic 
peak set limits on the electron scattering optical depth of a reionized IGM which 
imply z rc \ < 40 (model-dependent) [26]. Together with the lack of a hydrogen Lya 
resonance scattering ( "Gunn-Peterson" ) trough in the spectrum of an SDSS quasar 
discovered at z — 5.8 [19], these limits currently suggest that 6 < z Tci < 40. The 
origin and consequences of this reionization are among the major unsolved problems 
of cosmology. (For prior reviews and further references, see, e.g. [3,29,55].) Photons 
emitted by hitherto undetected massive stars or miniquasars formed within early 
galaxies are generally thought to be the reionization source. 

This is consistent with the theoretical expectation in cosmological models like 
the Cold Dark Matter (CDM) model, in which the first objects to condense out 
of the background and possibly begin star formation were small, of subgalactic 
mass (i.e. M ^ 10 6 M Q ) and began to form as early as z ^ 30 (e.g. [25,30,57,71]). 
The initial collapse of such objects only led to star formation if radiative cooling 
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was possible after collapse, usually involving H 2 molecules, and it is likely that the 
first objects to form stars released radiation (and possibly SN explosion energy as 
well) which exerted a strong feedback on the subsequent formation history of other 
objects. The details of this feedback and even the overall sign (i.e. negative or 
positive) are poorly understood (e.g. [8,14,17,20,28,31,32,36,55,57,59]). It appears 
that the first stars would have photodissociated H 2 long before enough UV was 
emitted to bring about reionization, so it is currently thought that if reionization 
was accomplished by stars, they formed in objects with virial temperature T vir > 
10 4 K, of mass > 10 8 M Q , which were able to cool by atomic radiative cooling, even 
without H 2 . This conclusion changes, however, if the first sources were miniquasars 
whose nonthermal spectra had a significant X-ray flux, since this would have created 
a positive feedback on the H 2 . 

If we adopt an optimum efficiency for massive star formation and radiation release 
by the collapsed baryon fraction in a standard CDM model which is flat, matter- 
dominated, and COBE-normalized (an unrealistic model which has too high an 
amplitude to satisfy X-ray cluster abundance constraints at z — 0, but conserva- 
tively overestimates z rei ), z Tci 50 is possible [55]. But more suitable CDM models 
all tend to yield z rci < 20 (e.g. [4,13,15,25,30,57,74]). 

II COSMOLOGICAL IONIZATION FRONTS 

Ionization Fronts in a Clumpy Universe. The neutral, opaque IGM out 
of which the first bound objects condensed was dramatically reheated and reion- 
ized at some time between a redshift z ~ 50 and z ~ 6 by the radiation released 
by some of these objects. When the first sources turned on, they ionized their 
surroundings by propagating weak, R-type ionization fronts which moved outward 
supersonically with respect to both the neutral gas ahead of and the ionized gas 
behind the front, racing ahead of the hydrodynamical response of the IGM [54,56]. 
The problem of the time-varying radius of a spherical I-front which surrounds iso- 
lated sources in a cosmologically-expanding IGM was solved analytically by [54,56], 
taking proper account of the I-front jump condition generalized to cosmological 
conditions. They applied these solutions to determine when and how fast these 
I-front-bounded spheres would grow to overlap and, thereby, complete the reion- 
ization of the universe. The effect of density inhomogeneity on the rate of I-front 
propagation was described by a mean "clumping factor" q > 1, which slowed the I- 
fronts by increasing the average recombination rate per H atom inside clumps. This 
suffices to describe the average rate of I-front propagation as long as the clumps 
are either not self-shielding or, if so, only absorb a fraction of the ionizing photons 
emitted by the central source. In two recent calculations, this analytical prescrip- 
tion was adapted to N-body simulations of structure formation in the CDM model, 
to calculate the evolving size of the spherical H II regions with which to surround 
putative sources of ionizing radiation, to model the growth of the ionized volume 
filling factor leading to reionization [4,15]. 



Numerical radiative transfer methods are currently under development to solve 
this problem in 3D for the inhomogeneous density distribution which arises as cos- 
mic structure forms, but so far are limited to an imposed density field without gas 
dynamics (e.g. [1,16,45,51]). A different approach, which replaces radiative transfer 
with a "local optical depth approximation," intended to mimic the average rate at 
which I-fronts expanded and overlapped during reionization within the context of 
cosmological gas dynamics simulation, has also been developed [22]. These recent 
attempts to model inhomogeneous reionization numerically are handicapped by 
their limited spatial resolution (> lkpc), which prevents them from resolving the 
most important density inhomogeneities. The dynamical response of density inho- 
mogeneities to the I-fronts which encountered them and the effect which these inho- 
mogeneities had, in turn, on the progress of universal reionization, therefore, require 
further analysis. Toward this end, we have developed a radiation-hydrodynamics 
code which incorporates radiative transfer and have focused our attention on prop- 
erly resolving this small-scale structure. In what follows, we summarize the results 
of new radiation-hydrodynamical simulations of what happens when a cosmological 
I-front overtakes a gravitationally-bound density inhomogeneity - a dwarf galaxy 
minihalo - during reionization. According to [27], the photoevaporation of these 
sub-kpc-sized objects is likely to be the dominant process by which ionizing pho- 
tons were absorbed during reionization, so this problem is of critical importance in 
determining how reionization proceeded. 

Dwarf Galaxy Minihalos at High Redshift. The effect which small-scale 
dumpiness had on reionization depended upon the sizes, densities, and spatial dis- 
tribution of the clumps overtaken by the I-fronts during reionization. For the 
currently-favored ACDM model (fi = 1 - A = 0.3, h = 0.7, Q b h 2 = 0.02, 
primordial power spectrum index n p — 1; COBE-normalized), the universe at 
z > 6 was already filled with dwarf galaxies capable of trapping a piece of the 
global, intergalactic I-fronts which reionized the universe and photoevaporating 
their gaseous baryons back into the IGM (see Figure 1). Prior to their en- 
counter with these I-fronts, "minihalos" with T V - 1T < 10 4 K were neutral and opti- 
cally thick to hydrogen ionizing radiation, as long as their total mass exceeded 
the Jeans mass Mj in the unperturbed background IGM prior to reionization [i.e. 
Mj = 5.7xl0 3 (fi /i 2 /0.15)- 1 / 2 (fi 6 /i 2 /0.02)- 3 / 5 ((l+2)/10) 3 / 2 M Q ],aswas required to 
enable baryons to collapse into the halo along with dark matter. Their "Stromgren 
numbers" L$ = 2R^ a i Q /is, the ratio of a halo's diameter to its Stromgren length is 
inside the halo (the length of a column of gas within which the unshielded arrival 
rate of ionizing photons just balances the total recombination rate), were large. For 
a uniform gas of H density nn,c, located a distance rM pc (in Mpc) from a UV source 
emmitting N p ^q ionizing photons (in units of 10 56 s _1 ), the Stromgren length is 
only £ s « (100pc)(iVp hi 5 6 /r 2 Jpc )(n ffiC /0.1 cm -3 )" 2 , so L s » 1 for a wide range of 
halo masses and sources of interest. In that case, the intergalactic, weak, R-type 
I-front which entered each minihalo during reionization would have decelerated to 
about twice the sound speed of the ionized gas before it could exit the other side, 
thereby transforming itself into a D-type front, preceded by a shock. Typically, the 



side facing the source would then have expelled a supersonic wind backwards toward 
the source, which shocked the surrounding IGM as the minihalo photoevaporated. 

The importance of this photoevaporation process has long been recog- 
nized in the study of interstellar clouds exposed to ionizing starlight (e.g. 
[5,6,38,39,41,43,47,53,69]). In the cosmological context, however, its importance 
has only recently been fully appreciated. 

In proposing the expanding minihalo model to explain Lya forest ( "LF" ) quasar 
absorption lines, [9] discussed how gas originally confined by the gravity of dark- 
matter minihalos in the CDM model would have been expelled by pressure forces if 
photoionization by ionizing background radiation suddenly heated all the gas to an 
isothermal condition at T ~ 10 4 K, a correct description only in the optically thin 
limit. The first discussion of the photoevaporation of a primordial density inho- 
mogeneity overtaken by a cosmological I- front, including radiation-hydrodynamical 
simulations, however, was in [65,66]. [2] subsequently estimated the relative impor- 
tance of this process for dwarf galaxy minihalos of different masses at different 
epochs in the CDM model, concluding that 50%-90% of the gas which had already 
collapsed into gravitationally bound objects when reionization occurred should have 
been photoevaporated. 

Not only did this photoevaporation during reionization affect most of the col- 
lapsed baryons, but so common were minihalos with T v ; r < lO 4 ^ during reioniza- 
tion that the typical reionizing photon is likely to have encountered one of them, 
according to [27]. A catalogue of all halos in a cubic volume 100 kpc (50 kpc) on a 
side at z = 9 in the currently-favored ACDM model based on N-body simulations 
reveals that minihalos with T vir < 10 A K (i.e. M < 1O 76 M ) were separated on 
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FIGURE 1. ACDM Halos at High Redshift: DM Halos at z = 9 are shown with sizes and 
locations determined by FOF algorithm applied to P 3 M simulations (128 3 particles, 256 3 cells), 
projected onto face of simulation cube of proper size Lbox, as labelled, (a) (left) Mhaio > 10 7 - 6 M Q 
only (i.e. T vir > 10 4 A). (b) (middle) 10 5 6 M Q < M halo < 10 7 6 M Q only (i.e. T vir < 10 4 A) - 
THE MINIHALOS; (c) (right) like (b), but higher resolution simulation with same number of 
particles and cells in 1/8 volume (i.e. M hal0imin = 10 4 - 7 M Q ). 



average by only d « 7kpc (3.5 kpc) for M > 1O 5 - 6 M (M > 1O 4 - 7 M ), respectively, 
while their geometric cross sections together covered / ~ 16% (30%) of the area 
along every 100 kpc of an average line of sight [60] [see Figs. 1(b), (c)]. If the 
sources of reionization, on the other hand, were larger-mass halos with T v - lT > 10 A K 
[like those in Fig. 1(a)], then these were well-enough separated that typical reion- 
ization photons were likely to have been absorbed by intervening photoevaporating 
minihalos. To demonstrate this in a statistically meaningful way with more dy- 
namic range than the N-body results can yet provide, an analytical approximation 
to the detailed numerical results is required. We shall combine the well-known 
Press-Schechter (PS) prescription for deriving the average number density of halos 
of different mass at each epoch with the nonsingular truncated isothermal sphere 
(TIS) model of [34,58] for the size, density profile, and temperature of each halo 
as unique functions of the halo mass and collapse redshift for a given background 
universe. (The TIS model and further tests and applications of it are briefly men- 
tioned elsewhere in this volume [35]). We illustrate the validity of these aproxima- 
tions in Figure 2 by comparing them with the N-body simulation results depicted 
in Figure 1. The PS halo mass functions plotted in Figure 2(a) are shown in Fig- 
ure 2(b) to reproduce the N-body results over a range of 10 5 in halo mass, for 
which dn halo /dM » lO 12 (M/M )" 2 Mpc _3 M (proper units) at z = 9. Likewise, 
Figure 2(c) shows that the TIS model correctly predicts the average virial ratio, 
[GM 2 oo/(c r y^2oo)]Tis = 2.176, for halos of different mass according to the N-body 
results, over this same mass range. 

This TIS+PS approximation allows us to determine which halo masses are sub- 
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FIGURE 2. CDM Halos: N-body Results Vs. PS+TIS Approximations: (a) (left) 
Proper Number Density of Halos per Unit Halo Mass, dn^io/dM , in ACDM at different rcdshifts, 
as labelled, according to PS approximation; (b) (middle) Same as (a), with PS mass function at 
z = 9 (solid curve) compared with N-body results for halos in Fig. 1 for three different resolutions: 
£box = 100 kpc (circles), 50 kpc (crosses), and 25 kpc (triangles); (c) Virial ratio, GM 2 oo/ (&v r 20o) , 
versus halo mass M 2 oo, where <jy is the halo DM velocity dispersion and M 2 oo and r 2 oo are the 
mass and radius of the sphere whose average density is 200 times the cosmic mean value, for all 
halos in Fig. 1 with at least 200 particles per halo [symbols same as in (b), with la error bars]. 
Horizontal line is analytical prediction of the TIS model. 




ject to photoevaporation and how common they £1X6 ; clS functions of redshift. Each 
minihalo with T v ; r < 10 4 K is opaque to H-ionizing photons with a geometric 
cross section (Xhaio — 7rrf , where r t = TIS radius ~ (0.75 kpc)(M/lO 7 M ) 1 / 3 [(l + 
-2coii)/10]~ 1 (fio^ 2 /0.15) -1 / 3 . Ionizing photons travelling through this universe will 
suffer absorption by these minihalos [i.e. those with mass Mj < M < M(10 4 i^) 
at each z] with a mean free path A m f p , as shown in Figure 3(a). For comparison, 
the mean separation < d scp > of halos of each mass is plotted in Figure 3(b). At 
z = 9, for example, A m f p = 160 kpc, while halos of mass M > 10 8 M Q are separated 
on average by <d scp >m 5O(M/lO 8 M ) 1//3 kpc. As shown in Figure 3(c), their ratio, 
< ^sc P > /A m fp, gives the fraction of the sky, as seen by a source halo of a given 
mass, which is covered by opaque minihalos located within the mean volume per 
source halo. If halos with M > 10 8 M Q are the reionization sources, their minihalo 
covering fraction is close to unity and increases with increasing source mass. This 
estimate will increase by a factor of a few if we take account of the statistical bias 
by which minihalos tend to cluster around the source halos [27]. 

An argument like this led [27] to argue that our photoevaporating minihalos 
are the chief consumers of the ionizing photons responsible for reionization. As a 
result, they suggest, the photoevaporation of these minihalos drives the number 
of ionizing photons per baryon required to reionize the universe up by an order 
of magnitude compared to previous estimates! A recent semi-analytical study of 
inhomogeneous reionization by [44], which neglected this effect, concluded that only 
one ionizing photon per hydrogen atom would have been sufficient to reionize most 
of the volume of the IGM by z ~ 5, an order of magnitude too low when compared 
to the new photoevaporation-dominated estimates. The gas clumping model they 
adopted apparently missed the smallest scales because it was adjusted to match 
numerical simulation results which could not resolve these scales. Simulations by 
[22], which agreed with [44], have the same problem since their resolution limit 




FIGURE 3. Minihalo Sinks and Sources: (a) (left) Proper mean free path <ni aa ioO'haio > _1 
for absorption of photons by minihalos in ACDM at different z (i.e. if they photoevaporate at 
z cv < z); (b) (middle) Proper mean separation of halos <d sop >= (M dn^aio/dM) -1 / 3 versus halo 
mass at different redshifts in ACDM, as labelled; (c) (right) Fraction of sky covered by minihalos 
located within the mean volume per source halo, F COVOINSOUrco =<d sop > /A m f p , versus source halo 
mass, at different redshifts. 



exceeded ~ 1 kpc, too large to resolve the minihalos which photoevaporate. 

We are led to conclude that further study of the photoevaporation of cosmolog- 
ical minihalos during reionization is essential if we are to advance the theory of 
reionization. The new "photoevaporation-dominated" reionization scenario sug- 
gests that significantly more than one photon per baryon may have been required, 
and this is difficult to understand on the basis of simple extrapolation (e.g. [21,42]) 
of either quasar or stellar photon production rates observed at z < 5 to z > 5. 
Perhaps this problem will be alleviated by appeal to the recently revised stellar 
output for zero metallicity stars [10,73], or to a higher escape fraction of ioniz- 
ing photons from their source than the nominal / csc ~ 0.1 typically assumed (e.g. 
[18,40,52,70,76]). Alternatively, if minihalo sources, alone, reionized the IGM (e.g. 
with miniquasars instead of starlight), then their ionizing photons could have done 
so without encountering any other minihalos, since -F C over,source < 1 for halos with 
T vir < 10 4 K. Photons released, thereafter, would have encountered fewer opaque 
minihalos, as well, since baryons would not then have condensed out of a reheated 
IGM to form new minihalos, according to [57]. 

The good news is that observations should be able to distinguish these possibili- 
ties. In order for photoevaporating minihalos to have consumed most of the ionizing 
photons produced before reionization is complete, the covering fraction of the ion- 
ization sources by these photoevaporative flows must be of order unity. Hence, if 
we can observe the sources of universal reionization directly, we will generally be 
able to observe foreground photoevaporating minihalos in absorption toward these 
sources. Such observations will constrain and help diagnose the reionization pro- 
cess. As shown in Figure 4, photons emitted by any high z source before or during 
reionization will typically encounter large numbers of photoevaporating minihalos 
at z > 6. The number of minihalos probed per unit redshift interval, in fact, is just 
(1 + z)^ 1 times the ratio of the horizon size c/H(t) at z to A m f p , a large number. 
It is important, therefore, for us to predict the effect which reionization will have 




FIGURE 4. Minihalos Encountered By a Photon Emitted At High Redshift: (a) (left) 
Number of minihalos per unit redshift interval in ACDM encountered by a photon at redshift z 
which travels along the average LOS; (b) (middle) Total number of minihalos in ACDM along 
photon path which probes redshift interval z m i n < z < z max for z m i n = 6 is plotted versus z max ; 
(c) (right) Asymptotic number of halos along photon path which probes redshift interval (z, oo). 



on the gas in these minihalos. 

The Photoevaporation of Dwarf Galaxy Minihalos Overtaken by Cos- 
mological Ionization Fronts. We have performed radiation-hydrodynamical 
simulations of the photoevaporation of a cosmological minihalo overrun by a weak, 
R-type I-front in the surrounding IGM, created by an external source of ionizing ra- 
diation [61-64]. The gas contained H, He, and a possible heavy element abundance 
of 1CT 3 times solar. Our simulations in 2D, axisymmetry used an Eulerian hydro 
code with Adaptive Mesh Refinement and the Van Leer flux-splitting algorithm, 
which solved nonequilibrium ionization rate equations (for H, He, C, N, O, Ne, and 
S) and included an explicit treatment of radiative transfer by taking into account 
the bound-free opacity of H and He [43,49,50]. The reader is referred to [65-67] for 
earlier results which considered uniform clouds and demonstrated the importance 
of a proper treatment of optical depth. 

Here we compare some of our results of minihalo photoevaporation in a ACDM 
universe [64] for two types of source spectra: a quasar-like source with emission 
spectrum F v oc z/~ L8 [y > v H ) and a stellar source with a 50,000 K blackbody 
spectrum, with luminosity and distance adjusted to keep the ionizing photon fluxes 
the same in the two cases. In particular, if ru pc is the proper distance (in Mpc) be- 
tween source and minihalo and N p ^ 56 is the H-ionizing photon luminosity (in units 
of 10 56 s _1 ), then the flux at the location of the minihalo would, if unattenuated, 
correspond initially to Np^^/r^^ = 1; thereafter, rM pc oc a(t), the cosmic scale 
factor. 

Our initial condition before ionization is that of a 1O 7 M minihalo in the ACDM 




FIGURE 5. Photoevaporating Minihalo I. One time-slice, 60 Myr after I-front caused by 
ionizing source (located far to the left of computational box along the a;-axis) overtakes a 1O 7 M0 
minihalo [centered at (r,x) = (0, 1.06 x 10 22 cm)] at z = 9 in the ACDM universe, for two types 
of source spectra: (a) (left) STELLAR CASE and (b) (right) QUASAR CASE, (upper panels) 
isocontours of atomic density, logarithmically spaced, in (r,x)— plane of cylindrical coordinates; 
(lower panels) velocity arrows are plotted with length proportional to gas velocity. An arrow of 
length equal to the spacing between arrows has velocity 25kms _1 ; minimum velocities plotted 
are 3kms _1 . Solid line shows current extent of gas initially inside minihalo at z = 9. Dashed line 
is I-front (50% H-ionization contour). 



universe which collapsed out and virialized at z co \\ = 9, yielding a truncated, nonsin- 
gular isothermal sphere ("TIS") of radius r t = 0.75 kpc in hydrostatic equilibrium 
with T vir = 4000 K and dark-matter velocity dispersion ay = 5.2 km s" 1 . The TIS 
profile has a central density and an average density which are 18,000 and 130 times 
the mean background density, respectively, with core radius r = 7"King/3 ~ r t /30. 
This hydrostatic sphere is embedded in a self-similar, spherical, cosmological infall 
according to [7]. 

The results of our simulations on an (r, x)-grid with 256 x 512 cells are illustrated 
by Figures 5-9. The minihalo shields itself against ionizing photons, traps the 
R-type I-front which enters the halo, causing it to decelerate inside the halo to 
close to the sound speed of the ionized gas and transform itself into a D-type 
front, preceded by a shock. The side facing the source expels a supersonic wind 
backwards towards the source, which shocks the IGM outside the minihalo. The 
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FIGURE 6. Photoevaporating Minihalo II. Same time-slice as Fig. 5. (a) (left) STEL- 
LAR CASE and (b) (right) QUASAR CASE. From top to bottom: (i) isocontours of pressure, 
logarithmically spaced, in (r,x)— plane of cylindrical coordinates; (ii) pressure along the r = 
symmetry axis; (hi) temperature; (iv) H I fraction; (v) He I (solid) and He II (dashed) fractions; 
(vi) bound-free optical depth measured from x = along the x-axis, at the threshold ionization 
energies for H I (solid), He I (dashed), He II (dotted). Key features of the flow are indicated by the 
numbers which label them on the temperature plots: 1 = IGM shock; 2 = contact discontinuity 
which separates shocked halo wind (between 2 and 3) from swept-up IGM (between 1 and 2); 
3 = wind shock; between 3 and 4 = supersonic wind: 4 = I-front; 5 = boundary of gas initially 
inside minihalo at z — 9; 6 = shock in shadow region caused by compression of shadow gas by 
shock-heated gas outside shadow. 



wind grows more isotropic with time as the remaining neutral halo material is 
photoevaporated. Since this gas was initially bound to a dark halo with ay < 
lOkms -1 , photoevaporation proceeds unimpeded by gravity. Figures 5 and 6 show 
the structure of the photoevaporative flow 60 Myrs after the global I-front first 
overtakes the minihalo, with key features of the flow indicated by the labels on 
the temperature plot in Figure 6. Figure 7(a) shows how the neutral mass of 
the gas initially within the original hydrostatic sphere gradually declines as the 
minihalo photoevaporates, within t ev w 250 (100) Myrs for the stellar (quasar) 
cases, respectively. The gradual decay of the opaque cross section of the minihalo 
as seen by the source is illustrated by Figure 7(b). 

Some observational signatures of this process are shown in Figures 8 and 9. 
Figure 8 shows the spatial variation of the relative abundances of C, N, and O 
ions along the symmetry axis after 60 Myrs. While the quasar case shows the 
presence at 60 Myrs of low as well as high ionization stages for the metals, the 
softer spectrum of the stellar case yields less highly ionized gas on the ionized side 
of the I-front (e.g. mostly C III, N III, O III) and the neutral side as well (e.g. C II, 
N I, O I and II). The column densities of H I, He I and II, and C IV for minihalo 
gas of different velocities as seen along the symmetry axis at different times are 
shown in Figure 9. At early times, the minihalo gas resembles a weak Damped Lya 
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FIGURE 7. Evolution of Neutral Gas Content of Photoevaporating Minihalo. 

(a) (left) Fraction of mass Mj, the mass which is initially inside the minihalo when the in- 
tergalactic I-front overtakes it, which remains neutral versus time t (in units of the evaporation 
time t ev at which this M/Mi = 1CP 3 ), for a large range of cases with different assumed source 
spectra and fluxes and minihalo masses; (b) (right) fraction of minihalo initial geometric cross 
section irrf which is opaque to source photons that can ionize H I (top panel), He I (middle panel), 
or He II (bottom panel) versus time (in units of t cv ), for the case with stellar source shown in 
Figs. 5, 6 (t cv = 245 Myr). 



("DLA") absorber with small velocity width (> lOkms" 1 ) and jV H i ^ 10 20 cm- 2 , 
with a Lya-Forest("LF")-like red wing (velocity width > lOkms -1 ) with iV H i ^ 
10 16 cm -2 on the side moving toward the source, with a He I profile which mimics 
that of H I but with Nn e i/Nm ~ [Hc]/[H], and with a weak C IV feature with 
iVciv ~ 10 11 (10 12 )£ cm -2 for the stellar (quasar) cases, respectively, displaced in 




FIGURE 8. Observational diagnostics I: ionization structure of metals. C, N, and O 

ionic fractions along symmetry axis at t — 60Myr, for photocvaporating minihalo of Figs. 5, 6. 
(a) (left) STELLAR CASE; (b) (right) QUASAR CASE. 
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FIGURE 9. Observational diagnostics II. Absorption lines. Minihalo column densities 
(cm -2 ) along symmetry axis for gas at different velocities, for photoevaporating minihalo of 
Figs. 5, 6. (a) (left) STELLAR CASE; (b) (right) QUASAR CASE. (Top) H I; (Middle) He I 
(solid) and He II (dotted); (Bottom) C IV (i.e. if [C]/[C]© = £ x 10~ 3 , then plotted values are 
Nctv/0- Each box labelled with time (in Myrs) since arrival of intergalactic I-front. 



this same asymmetric way from the velocity of peak H I column density, where 
£ = [C]/[C]© x 10 3 . For He II at early times, the stellar case has iV He ii ~ 10 18 cm~ 2 
shifted by 10's of km/sec to the red of the H I peak, while for the quasar case, 
He II simply follows the H I profile, except that N^n/Nm ~ 10 in the red wing 
but Nueii/Nm ~ 10~ 2 in the central H I feature. After 260 (160) Myr, however, 
only a LF-like H I feature with column density N m ~ 10 14 cm~ 2 remains, with 
AW^hi ~ 1/4 (< 1/30), N Uell /N m ~ 10 3 (10 2 ), and AWiVui ~ 3(1) x[C]/[C] 
for the stellar (quasar) cases, respectively. 

As described above, intervening minihalos like these are expected to be ubiquitous 
along the line of sight to high redshift sources. With photoevaporation times t ev > 
100 Myr, this process can continue down to redshifts significantly below z = 10. For 
stellar sources in the ACDM model, these simulations show that photoevaporation 
of a 10 Mq minihalo which begins at ^initial = 9 can take 250 Myr to finish, at za n3i \ = 
6.8, during which time such minihalos can survive without merging into larger 
halos. Observations of the absorption spectra of high redshift sources like those 
which reionized the universe should reveal the presence of these photoevaporative 
flows and provide a useful diagnostic of the reionization process. 

Dwarf Galaxy Suppression and Reionization. As emphasized by [57], the 
reheating of the IGM which accompanied its reionization must have filtered the 
linear growth of baryonic fluctuations in the IGM, thereby reducing the baryon 
collapsed fraction and preventing baryons from condensing out further into the 
minihalos (see also [12,13,23,24,74]). A related effect, the suppression of baryon 
accretion onto dark matter halos, has also been studied, by ID [33,36,37,72] and 
3D simulations (e.g. [46,48,75]). The current conclusion is expressed in terms of the 
threshold circular velocity, v c , below which baryonic infall and star formation were 
suppressed by photoionization: v c ~ 30kms _1 for complete suppression, partial 
suppression of infall found to extend even to v c ~ 75kms~ 1 . This suppression 
of gas accretion onto low-mass halos by the feedback effect of reionization may 
naturally explain why there are so many fewer dwarf galaxies observed in the Local 
Group than are predicted by N-body simulations of the CDM model [11]. 
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